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ABSTRACT 

We report on TeV 7-ray observations of the blazar Mrk421 (redshift of 0.031) with the 
VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of 
VERITAS allowed us to sample the TeV 7-ray fluxes and energy spectra with unprecedented 
accuracy where Mrk421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS 
and 96 hrs of Whipple 10 m data were acquired between January 2006 and June 2008. We present 
the results of a study of the TeV 7-ray energy spectra as a function of time, and for different flux 
levels. On May 2nd and 3rd, 2008, bright TeV 7-ray flares were detected with fluxes reaching 
the level of 10 Crab. The TeV 7-ray data were complemented with radio, optical, and X- 
ray observations, with flux variability found in all bands except for the radio waveband. The 
combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux 
levels, often correlated with an increase of the source activity in TeV 7-rays. Contemporaneous 
spectral energy distributions were generated for 18 nights, each of which are reasonably described 
by a one-zone SSC model. 

Subject headings: BL Lacertae objects: individual (Mrk421) — Galaxies: jets — Galaxies: nuclei — 
Gamma rays: galaxies — X-rays: galaxies 
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1. Introduction 

In 1992, observations with the Whipple 10m 
Cherenkov telescope led to the first discovery of 
an extrag alactic source of Te V 7-rays, the blazar 
Mrk421 (jPunch et al. I Il992h - Since then, more 



than 30 similar sources have been detected with 



groun d-based 7-ray detectors (jWakelv fc Horan 
201df ). The sources with well measured red shifts 



lie from 0.031 (Mrk 421. [Snellen et al. I ([20021) ) to 



0.536 for the recently detected radio quasar 3C 279 
([Albert et al. I 120081) . Typically, blazars show 
core-dominated emission, and they are character- 
ized by rapid variability. Their spectral energy dis- 
tribution (SED) in the vF v representation is char- 
acterized by two broad, well-separated "humps" 
arising from (i) synchrotron emission (low-energy) 
and (ii) a high-energy component of either lep- 
tonic or hadronic nature. During TeV 7-ray 
flares, strong sources (i.e. Mrk 421 and PKS2155- 
304) exhibit ^i^-fluxes of about 10~ 9 erg cm" 2 s _1 



([Aharonian et al. I l2009bl) ; corresponding 



to 7- 

ray luminosities of between 10 42 and lO^ergs" 1 
for assumed anisotropic emission with an open- 
ing angle of 5 deg. The blazars detected at TeV 
energies are the high-frequency peaked counter- 
parts of the blazar source population detected at 
MeV/GeV energies with the EGRET experiment 
on board the space-borne Compton Gamma Ray 
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Observatory ( Hartman et al. 1999), and recently 
expanded by the Large Area Tele scope (LAT) 
on b oard the Fermi 7-ray satellite ( Abdo et al. I 
20101) . 

Following the detection in TeV 7-rays in 1992, 
Mrk421 was observed intensively, and the obser- 
vations led to a number ol landmark discoveries: 

1. The detection of fast flux v ariability with 



a do ubling time of 15min (jGaidos et al 



1996). 



2. The first tentative evi dence for a X-ray/TeV 



7-ray flux correlation ( Buckley et al. 19961 ) 



Several observation campaigns strengthened 
the evidence for such a correlation. Some of 
th e most convinc i ng re sults were presented 



bv lFossati et al. I (2008). 



3. The detection of TeV 7-ray energy spec- 
tra that harden with increasing fluxes 
(|Krennrich et al. II2OO2} ). 



Modeling of Mrk421 data with synchrotron- 
Compton models revealed the first evidence for 
bulk jet Lorentz factors of the order of 50 



( Krawczvnski et al~ll200ll) and modeling of data 
taken during different states revealed evidence 
that one-zone synchrotron self-Compton (SSC) 
model s are insufficient to des cribe the observa- 
tions (jBlazeiowski et al. 2005). Alternative mod- 



els are discu ssed in the literature. Rec e nt pa - 
pers include Katarzynski k Walczewska I ( 2010h : 

(2008); 



Bnu.clK- r k Permed d2010h 
^pMiffvl (l2009h : 



Gao et al 



Lichti et al. 



Tammi 

Stecker et al.1 (|2007t ). These models involve more 
complicated geometries which lead to a larger 
number of free model parameters than in the SSC 
model. 

Mrk421 has been a frequent target of multi- 
wavelength campaigns. During some of its very 
short flares t he X-ray and TeV 7-r ay fluxes tracked 



each other ( Fossati et al. 1 120081 ). However, X 



ray flares that are not accompanied by TeV 7- 
ray flares and vice versa have also been o bserved 



(|Rebillot et al~ll2006l: iFossati et al. II2008D . There 
is good evidence that the X-ray and TeV 7- 
ray activities are correl ated when averaged over 



~1 week time inter vals (jBlazeiowski et al. 1 120051: 
Horan et al. 20091 ). All attempts to establish 



convincing evidence for a correlation of the X- 
ray and 7-ray fluxes with the flux variability at 



radio to optical wavelengths, have failed s o far 
(jBlazeiowski et al. II2OO5I : iHoran et al. Il2009h . 



The X-ray/TeV 7-ray correlation properties 
are of great interest, as they might enable us 
to decide on the emission model, i.e. between 
(i) leptonic models in which a single popula- 
tion of high energy electrons emits the low and 
high-energy radiation as synchrotron and inverse 
Compton (IC) emission, respectively, and (ii) 
hadronic models in which the 7-ray emission is 
attributed to additional population(s) of high- 
energy particles, powered by the acceleration of 
extremely high-energy protons. It is interesting 
to note that the X-ray/TeV 7-ray flux correla- 
tion has only been studied for a handful of sources 
(Mrk421, Mrk501, 1ES1959+650, PKS 2155-305, 
1ES2344+514) with a sufficiently good signal-to- 
noise ratio in both bands to investigate the cor- 
relation properties. Although the X-ray and TeV 
7-ray fluxes seem to be correlated, it is not clear 
how well this correlation holds for individual flares 



([Krawczvnski et al. II2004I ). 



In this paper we report on the first TeV 7-ray 
observation campaign on Mrk421 performed with 
the VERITAS observatory. VERITAS achieves 
an energy flux sensitivity of (£ o bs/50hrs) x 
5 x lO" 13 ergs cm" 2 sec" 1 at 500 GeV Q. The 
high sensitivity of VERITAS and the brightness 
of Mrk421 during flares allow us to measure fluxes 
with ^minute time bins, and to determine energy 
spectra for 5-min time intervals or less. We also re- 
port accompanying observations in the radio band, 
at optical wavelengths, and in the X-ray band with 
the Rossi X-ray Timing Explorer (RXTE) PCA, 
Swift XRT and Suzaku instruments. 

In Section [5] we discuss the data sets and the 
analysis methods. The results of different studies 
are presented in Section [3] followed by a summary 
and discussion in Section |U 

2. Data Sets and Data Reduction 

In this section we describe the observations 
and analysis of the data taken in the TeV en- 
ergy regime (Section [2. 1|) . at X-ray energies (Sec- 
tion l2.2|) . as well as in the optical (Section l2.3[) and 



iln terms of the flux from the Crab Nebula, the VERITAS 
sensitivity is (t o t, a /30min) -1//2 X 8% Crab in 30min of 
observations. This number is valid for the VERITAS array 
as operated in 2008. 
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radio ( Section 12.41) wavebands. 

2.1. VERITAS /Whipple 7 -ray data 

VERITAS. VERITAS consists of four 12 m di- 
ameter imaging atmospheric Cherenkov telescopes 
(IACTs) and is located at the base camp of the 
Fred Lawrence Whipple Observatory (FLWO) in 
southern Arizona at an altitude of 1280 m. It de- 
tects the Cherenkov light emitted by an extensive 
air shower (initiated by a 7-ray photon or cosmic 
ray entering the Earth's atmosphere) using a 499- 
pixel photomultiplier camera located in the focal 
plane of each telescope. The array is sensitive 
to 7-rays in the energy range from ^100 GeV to 
~30 TeV. Observations are performed (in 20min 
runs) on moonless nights using the "wobble" mode 
of operation, where all telescopes are pointed to 
a sky position offset of ±0.5 deg (alternating in 
direction between consecutive data runs) with re- 
spect to the source position. This method allows 
for a simultaneous background estimation to be 
made. More details about VERITAS, the data cal- 
ib ration and the analys is techniques can be found 



Acciari et al. I (|2008l) 



Only shower images which pass certain quality 
cuts are considered in the event reconstruction: 
image size > 500 digital counttH (dc) and image 
distance to the center of the camera <1.43 deg. 
The standard cuts for 7/hadron separation, which 
are base d on the width and l ength of the recorded 
images ( Acciari et al. |[2008h . were a priori opti- 
mized on data taken from the Crab nebula. An 
event is considered to fall into the signal (ON) re- 
gion once the squared angular distance between 
the reconstructed event direction and the Mrk421 
position is Ad 2 < 0.025 deg 2 . The background 
is estimated from different regions of the same 
size positioned at the same radial distance to 
the camera center as the ON region, and is re- 
ferred to as the refl ected background region model 
Berge et al. I (|2007l ). The excess is then calculated 
as the number of ON source counts less the nor- 
malized number of OFF region counts. In this 
analysis 5 OFF regions were used. The statisti- 
cal significances are calculated using the method 



2 Very Energetic Radiation Imaging Telescope Array System. 
3 The photomultiplier pulses are integrated within a time 

window of 24 ns duration. One photoelectron corresponds 

to approximately 5 digital counts. 



of iLi fc Ma I (Il983h . 

The energy Ei of an individual image is esti- 
mated using look-up tables generated from Monte 
Carlo simulations of 7-ray air showers. The ta- 
bles are parameterized in (i) the integrated charge 
(size), (ii) the impact parameter pi between the 
reconstructed shower axis and the optical axis of 
telescope i, (iii) the zenith angle z, (iv) the az- 
imuth angle Az, and (v) the level of the night-sky 
background. The energy of the shower event is 
then averaged over the n telescope energies to ob- 
tain -E r eco = l/ 71 Sr=i Ei- The energy resolution 
is estimated based on Monte Carlo simulations to 
be AE/E ?» 20% for energies between 100 GeV 
and 30 TeV. 

The effective areas A e g describe the energy- 
dependent response of the detector and are also 
obtained by Monte Carlo simulations. The effec- 
tive area A e g (E ieco , z, Az, AR, iVt e i) is estimated 
for each event based on its corresponding param- 
eters; AR is the angular distance between the 
reconstructed shower direction and the telescope 
pointing position, N tc i is the number of telescopes 
in the systerrQ. 

The inverse effective areas are used on an event- 
by-event basis to calculate the differential photon 
flux for each bin j in the energy spectrum: 



dN 
dE~ 



N, 



Tf vc ■ AEj \ f-j AeSM 



off J 



~J Aiff.k 
k—1 



(i) 

The sum of N on corresponds to reconstructed 
events from the ON region and the sum of N Q g to 
events from the five reflected OFF regions with the 
normalization a = 0.2. Tj lvc is the live time for bin 
j. The bias energy describes the energy at which 
the reconstructed energy E xeco (on average) devi- 
ates less than a certain percentage from the true 
energy Emc- The energy bias is calculated based 
on Monte Carlo simulations; a bias threshold of 
(E leco — Emc) /Emc = 10% is used in this analy- 
sis. The bias energy depends on the zenith angle 
of the corresponding data run. For each 20 min 
run, the bias energy is calculated and only those 
bins, j, in the energy spectrum which are fully 



4 A part of the Mrk421 data were taken with only N tc \ = 3 
operating telescopes. 
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contained above the bias energy are allowed to re- 
ceive events. The live time Tj lve is increased on a 
run- by-run basis for only those energy bins. 

In order to account for spill-over effects, the 
effective areas are calculated using the recon- 
structed energies, where the Monte Carlo input 
spectrum is weighted according to the recon- 
structed/measured spectrum in an iterative pro- 
cedure. The systematic errors of the parameters 
describing a power-law energy spectrum dN/dE = 
Io • E~ r have been estimated_| based on a Crab-like 
energy spectrum (r w 2.5) to be: AT/T = 8% and 
AJo/Io = 20%. 

The integral flux on a run-by-run basis above a 
certain energy E 1 (as shown in the light curves) is 
calculated as follows: A spectral slope is assumed 
and the effective areas for the corresponding run 
parameters (zenith angle, etc.) are used together 
with the measured excess to determine the nor- 
malization. The normalization is then used to cal- 
culate the integral flux above E' . This procedure 
has the advantage that the full event statistics are 
used and E' is not limited by the strongly vary- 
ing thresholds of individual data subsets (i.e. runs 
taken at different zenith angles). However, a spec- 
tral shape has to be assumed, which in our case 
is chosen (iteratively) for each data point accord- 
ing to the energy spectrum corresponding to the 
estimated flux level. 

Mrk421 is one of the objects in a trigger agree- 
ment between the 7-ray observatories H.E.S.S., 
MAGIC, and VERITAS put into place for known 
TeV 7-ray blazars in order to exchange informa- 
tion about flaring sources. The trigger criterion for 
Mrk421 is defined by a flux level measured by one 
of the observatories exceeding a value of 2$ cra b; 
this criterion was met several times during the 
campaign in 2008, leading to triggers sent by VER- 
ITAS and to triggers received from the MAGIC 
collaboration. VERITAS observed Mrk 421 during 
January/February and November/December 2007 
(5h) as well as in January- June 2008 (42.3 h) for 
a total of 47.3 h after run quality selection. The 
observation time corrected for the detector dead 
time amounts to 43.6 h. The zenith-angle range of 
the observations was 6deg— 56deg with an aver- 



age of 23.5 deg, corresponding to an analysis en- 
ergy threshold^ of 260 GeV. 

Whipple. The 10 m 7-ray Telesco pe at the Fred 



Lawr ence Whipple Observatory (jKildea et al 
20071 ) is sensitive in the energy range from 
200 GeV to 20 TeV with a peak response en- 
ergy (for a Crab-like spectrum) of approximately 
400 GeV. This telescope, although a factor of 
seven less sensitive than VERITAS, was used in 
this program to extend the TeV coverage when 
VERITAS was not available for Mrk 421 obser- 
vations. More detailed descriptions of Whipple 
observing modes and analysis procedures can be 
found elsewhere (Weekes 19961 : Punch k. Fegan 



199lt iRevnolds et al. I Il993l) . Details about the 



Whipple telescope includi ng the GRANITE- III 



camera have been given in iKildea et al. I (|2007j ) . 



The Whipple observations were conducted be- 
tween November 2005 and May 2008 (MJD 54417- 
54622). Only runs which pass the run quality se- 
lection (stability of the raw trigger rate, induced 
by the cosmic ray background) are considered in 
the analysis resulting in a data set of 96 h of ON- 
source data. The data were analyzed using the 
standard 2 nd - moment-parameterization technique 
(|Hillasl[l985h . Standard cuts (SuperCuts2000) 



were used to select 7-ray event s and to suppress 
backg round cosmic-ray events ( de la Calle et al. 1 
20031 ). Using the zenith angle dependence of a 
Crab data set taken at similar epochs, we account 
for the zenith angle dependence of the 7-ray excess 
rate by normalizing our measured Mrk 421 rate 
to the Crab rate at a corresponding zenith angle. 
It should be noted that this simplistic scaling is 
strictly only valid for a TeV spectrum close to that 
of the Crab Nebula (spectral index of V — 2.5). 
However, the systematic error introduced by this 
scaling can be expected to be small compared to 
the statistical error of the flux points. 

2.2. RXTE, SWIFT and Suzaku X-ray Ob- 
servations 

X-ray data were taken with the telesc opes on 
the R ossi X-r ay Timing Explorer R XTE \ Swank 
1994) , Swift (|Gehrels et al. 1 12004 ), and Suzaku 



5 Uncertainties in the atmosphere, components of the detec- 
tor, and shower reconstruction algorithms were considered 
in this estimate. 



6 The energy threshold is defined as the energy corresponding 
to the peak detection rate for a Crab-like spectrum. 
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(jMitsuda et al. 1 120071 ) satellites. The pointed X- 
ray observations are summarized in Table [TJ 

RXTE /PC A. The Propo rtional Counter Array 
(PCA) (jjahoda et al. Ill996l) comprises 5 Propor- 
tional Counter Units (PCUs) covering a nominal 
energy range of 2— 60 keV with a net detection area 
of 6250 cm 2 . Data between the energies of 3—15 
keV were used in this analysis. The 15 — 250 keV 
data from the H igh- Energy X-ray Timi ng Exper- 



Swift/XRT. The X-r ay telescope (XRT ) on 



board the Swift satellite (jGehrels et al. Il2004j ) is a 



iment HEXTE ( Rothschild et al. 1 119981 ) were not 



used owing to an insufficient signal-to-noise ra- 
tio. The PCA data were taken as part of a 
multi-wavelength observation proposal and com- 
prise 161 exposures between January 2006 and 
May 2008 (see Table P with a total net expo- 
sure time of 245.6 ksec. The observations had a 
typical exposure of 10 — 70 min per pointing and 
were taken at near-simultaneous times to sched- 
uled VERITAS observation^. For the observa- 
tions from January 6, 2006 to April 18, 2006 both 
PCU0 and PCU2 detectors collected data, while 
for all other data only PCU2 was operational. 
The data were filtered following the standard cri- 
teria advised by the NASA Guest Observer Facil- 
ity (GOF$. Standard-2 mode PCA data gathered 
with the top layer (X1L and X1R) of the opera- 
tional PCUs were analyzed using the HEAsoft 6.4 
package. Background models were generated with 
the tool pcarsp, based on the RXTE GOF calibra- 
tion files for a "bright" source with more than 40 
counts/sec. Response matrices for the PCA data 
were created with the script pcarsp. The saextrct 
tool was used to extract all PCA energy spectra. 

RXTE /ASM. The RXT E All Sky Monitor 
(ASM) (jLevine et al. 19961 ) is sensitive to X-ray 
energies at 2 — 12 keV and scans most of the sky 
every 1.5 hours. The data were obtained from the 
public MIT archiveEl in the form of 1-day aver- 
aged binning, as well as the dwell-by-dwell binning 
(for the short-term light curve and flux correlation 
studies). 



7 The RXTE/PCA staff at NASA GSFC and Principal Inves- 
tigator (PI) of the observation proposal Henric Krawczyn- 
ski, together with the VERITAS team, coordinated the ob- 
servations. 

8 http: / /heasarc. gsfc.nasa.gov/docs/xte/xte_ 1st. html 
9 http: / /heasarc. gsfc.nasa.gov/docs/xte/asm_ products. html 



focusing X-ray telescope w ith a 110 cm 2 effectiv e 



area and a 23' field of view ( Burrows et al. 2005) 



It is sensitive to X-rays in the 0.2 — lOkeV band. 
A total of 1175.7 ksec of XRT data were taken be- 
tween January 2006 and May 2008 (Table [jj in 
the Windowed Timing (WT) mode with grades 0- 
2 (referring to the pattern of CCD pixels for each 
event) selected over the energy range 0.4 — 10 keV. 
The XRTPIPELINE tool was used to calibrate 
and clean all Swift XRT event files with current 
calibration files. The data were reduced using 
the HEAsoft 6.4 package. Source counts were ex- 
tracted from a rectangular region of 40 pixels (94.4 
arcsec) along the one dimensional stream, and 20 
pixels high centered on the source. Background 
counts were extracted from a nearby source-free 
rectangular region of equivalent size. Ancillary re- 
sponse files were generated using the xrtmkarftask 
applying corrections for the PSF losses and CCD 
defects. The latest response matrix from the XRT 
calibration files was used. The extracted XRT en- 
ergy spectra were re-binned to contain a minimum 
of 20 counts in each bin and were fit with XSPEC 
12.4. 

Swift/BAT. The Burst Alert Telescope (BAT) 
is a large field of view (1.4 steradians) X-ray 
telescope with imaging cap abilities in the energ y 



range from 15 — 150 keV (jGehrels et al. 1 1200 



The BAT typically observes 50% to 80% of the 
sky each day. The data are the Swift/BAT tran- 
sient monitor results p rovided by the Swift/BAT 
team (jKrimm I l2008al) . Full details of the BAT 
data analysis are give n at the BAT transients web 
page (jKrimm I l2008bh . 



Suzaku/XIS . The X-ray Imaging Spectrometer 
(XIS) (jKovama et al. 1120071) on board the Suzaku 
satellite is composed of 3 X-ray CCD cameras 
combined with a single X-ray Telescope (XRT) 
covering a nominal energy range of 0.5 — 12keV. 
Each CCD camera covers an 18' x 18' region of the 
sky. The XIS data include observations between 
May 5, 2008 and May 8, 2008 with an exposure 
time of 180.8 ksec. Standard data reduction and 
processing were performed using HEAsoft v6.6.3 
and ftools v6.6. XIS events were extracted from 
a source region with an inner radius of 35 pix- 
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Table 1: Pointed X-ray observations 
(RXTE/PCA, Swift /XRT, and Suzaku/XIS) 
of Mrk421 in 2006 to 2008. 



curves whereby a unique flux offset is found for 
each spectral band (R, B, V) of eve ry instrument 



Start 


Stop 


N oha 


ObsID 


RXTE/PCA 


2006-01-06 


2006-03-02 


27 


91440-01 


2006-03-03 


2006-05-31 


48 


92402-01 


2008-01-07 


2008-05-07 


86 


93133-02 


Swift /XRT 








2006-01-02 


2006-12-05 


20 




2007-03-23 


2007-12-31 


24 




2008-01-07 


2008-05-08 


54 




Suzaku/XIS 








2008-05-05 


2008-05-08 


1 





els and an outer radius of 408 pixels. The ex- 
tent of the inner radius is such that pile-up ef- 
fects were minimized for the selected events. The 
background was selected from an annulus outside 
of the source region defined by 432 pixel and 464 
pixel inner and outer radii respectively. The re- 
sponse matrix and effective area were calculated 
for each XIS sensor using th e Suzaku ftools tasks , 
xisrmfgen and xissimarfgen (jlshisaki et al. 2007). 
XIS1 data were not included in this analysis. As 
the XIS0 and XIS3 have similar responses, their 
data were summed. 



2.3. Optical Observations 

Many optical observatories contributed data 
sets to this campaign (see below). The data from 
the observatories were reduced and the photome- 
try performed independently by different analysts 
using different strategies. The same set of ref- 
erence stars was used for all optical data sets to 
calculate the systematic error on the flux. How- 
ever, combining the various optical data to pro- 
duce a composite light curve for each spectral 
band is complicated by the fact that different 
observatories use different photometric systems. 
Furthermore, photometric apertures and the def- 
inition of the reported measurement error for 
each nightly-averaged flux is inconsistent across 
datasets. Therefore we have adopted a simple ap- 
proach for the construction of the composite light 



based on overlapping observations ([Steele et al. 
2007h . and the light curves have been scaled ac- 
cordingly (in our case the light curves of the Brad- 
ford Robotic Telescope and the New Mexico Skies 
observatory by 15% each). 

UVOT. Mrk421 was observed with the Swift 
Ultraviolet/Optical Telescope (UVOT) during 
2008. The instrument cycled through each of 
three ultraviolet pass bands, UVW1, UVM2 and 
UVW2 with central wavelengths of 260nm, 220nm 
and 193nm, respectively. More than 100 observa- 
tions were obtained with a typical/average expo- 
sure (per filter) of 150 s, ranging from 50 s up to 
900 s. Data were taken in the image mode, where 
the image is accumulated on board the satellite 
discarding the photon timing information within 
each single exposure to reduce the telemetry vol- 
ume and the time of transmission. 

Primary and secondary analyses were carried 
out using UVOTSOURCE standard tool and a 
custom UVOT pipeline. Both analyses used the 
calibration database released on February 2010. 
Photometry was computed using a 5 " source re- 
gion around the source and photometric c orrec 



tion s were ap p lied f ollowing iPoole et al. 
and Li et al. I ( 20061 ). All observations were in- 



(2008) 



spected manually. Astrometric misalignment be- 
tween the observed position and the nominal po- 
sition of Mrk421 which were found in several data 
sets were corrected by using a spatial fitting al- 
gorithm. Results of the two analysis chains were 
found to be in agreement. 

Due to the ultraviolet spectrum, we adopted 
the rate-to-flux conversion factors for GRB- 
like objects and not the standard factors 
used for Pickles-like star spectra. The fluxes 
were correc ted for galactic extin ction Eb-v — 



To obtain 
nearer values of 



0.015 mag (|Schlegel et al. I 119981 ) 
this value we smoothed the 
the database in correspondence with the coordi- 
nates of the source. Then, the computed opti- 
cal/U V galactic extinction coeffic ients were ap- 



plied (|Fitzpatrick fc Massa I Il999t l. The effects 



of intergalactic absorption and the zodiacal light 
have been estimated to be negligible and were not 
corrected for in this analysis. The fluxes and corre- 
sponding frequencies shown in the light curves are 
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redshift corrected, including a second-order cor- 
rection taking into account filter non-linearities. 

The host galaxy correction was not applied, 
but a systematic error on the flux is estim ated 
The measurements of Nilsson et al 



(1999) 



are 



used to estimate the host galaxy emission in the 
R band. These are used to obtain the corre- 
sponding components for the V, B and U bands 



(jFukugita et al. 1 11995). The presence of an up- 



turn flux excess in the far UV spectrum of ellipti- 
cal galaxies is caused by an old population of hot 
helium-burning stars without extended hydrogen- 
rich envelopes (as c ompared t o rat her young 
stars). The findings of lArimo"to" (|l996l ) were used 
to calculate the metallicity of the Mrk421 host 
galaxy and therefore constrain its con tribution to 
the UV bands to be less than 5% (|Han et al 



2007) 



Some caveats have to be mentioned. There are 
bright sources in the field of view which will cause 
significant coincidence losses, ghosting from inter- 
nal reflections, and may lead to an additional over- 
estimation of the Mrk421 blazar flux. Although 
the relative photometry (light curves) is expected 
to be less sensitive to these effects, an additional 
systematic error of 20% was added to the absolute 
UVOT fluxes shown in the SEDs to account for 
these uncertainties. 

BRT/NMS. Optical data were taken with the 
Bradford Robotic Telescope (BRT) in Tenerife, 
Canary Islands, Spain, as well as the New Mexico 
Skies observatory (NMS). The data were reduced 
by standard aperture photometr\F°l. The aperture 
size used was 10" diameter, and the comparison 
stars were taken from Villata et al. (1998). 

Bell. The Western Kentucky University's Bell 
observatory is a 0.6 m telescope located 12 miles 
southwest of Bowling Green, Kentucky. The ob- 
servations presented here were obtained with an 
AP6 CCD camera and Bessell R band filter. Dark 
and flat field corrections were made to the im- 
ages and differential aperture CCD photometry 
was performed using stars 1,3,2 of the compari- 
son sequence from Villata et al. (1998). No cor- 
rection for host galaxy flux or galactic absorption 



was made. 

WIYN. The WIYN 0.9m telescope is located 
at the National Optical Astronomy Observatory 
at Kitt Peak and operated by a consortium of 
universities. Observations of Mrk421 were per- 
formed since January 2006 in Johnson B and V 
and Cousins R optical filters and a CCD 1 deg field 
of view Mosaic Imager. Image reduction was per- 
formed with IRAF, using bias frames and dome 
flat fields for each night of data. We obtained 
magnitudes by differential photometry, using three 
reference stars from Villata et al. (1998). Since we 
are mostly interested in measuring the magnitude 
relative variations with time, these data were not 
corrected for the host galaxy flux or absorption. 

Tuorla/KVA. The Kungliga Vetenskapsaka- 
demien telescope (KVA, Royal Swedish Academy 
of Sciences) is located on Roque de los Mucha- 
chos, La Palma and operated by the Tuorla Ob- 
servatory, Finland. The telescope is composed 
of a 0.6 m f/15 Cassegrain devoted to polarime- 
try, and a 0.35 m f/11 SCT auxiliary telescope for 
multicolor photometry. This telescope has been 
successfully operated in a remote way since au- 
tumn 2003. Mrk421 has been observed in op- 
tical R-band typically once per night. Photo- 
metric measurements were made in differential 
mode, i.e. by obtaining CCD images of the target 
and calib rated comparison stars i n the same field 



of view dFiorucci fc Tosti 
19981: IVillata et al. Hl998l)T 



19961: iFiorucci et al. 



10 The MIRA Pro Version 7 (Mirametrics, Inc.) was used, see 
http: //www.mirametrics . com/ 



2.4. Radio Observations 

Radio data presented here were taken at four 
frequencies at two different radio observatories. 
The fluxes are given in Janskys (Jy), so they have 
already been normalized for the bandwidth of their 
receivers. 



Metsahovi The Metsahovi radio telescope 
(radome enclosed paraboloid antenna, diameter 
of 13.7 m) is situated in Finland. The measure- 
ments were made with a 1 GHz-band dual beam 
receiver centered at 36.8 GHz. The observations 
are ON-ON observations (typical integration time 
of 1200 — 1400 s), alternating the source and the 
sky in each feed horn. 
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The detection limit of the telescope at 37 GHz 
is of the order of 0.2 Jy under optimal conditions. 
Data points with a signal-to-noise ratio < 4 are 
handled as non-detections. The flux density scale 
is set by observations of DR 21. Sources 3C 84 and 
3C274 are used as secondary calibrators. A de- 
tailed descrip tion of the data reduction and analy- 
sis is given in Terasranta et al.1 ( 1998 ) . The error 



the flux and spectral variability in other energy 
bands in Section [3T4l 



estimate in the flux density includes the contri- 
bution from the measurement rms and the uncer- 
tainty of the absolute calibration. 

UMRAO. The University of Michigan Ra- 
dio Astronomy Observatory (UMRAO) (26 m 
paraboloid) provided monitoring data of Mrk421 
at 4.8 GHz, 8 GHz and 14.5 GHz between June 
2006 and May 2008. Each observation consisted 
of a series of ON-OFF measurements taken over a 
30-40 minute time period. All observations were 
made within a total hour angle range of about 5 
hours centered on the meridian. The calibration 
and reduction pro cedures have been described in 



3.1. Light Curves 

The radio, optical, X-ray and TeV light curves 
of Mrk421 are shown in Figure Q] for the years 
2006-2008 during which the source was extensively 
monitored in the various energy bands. A zoomed 
version for the year 2008 is shown in Figure[2j Two 
further levels of zoom are shown in Figure [3] (the 
two more active X-ray/TeV states) and Figured] 
(the strong TeV 7-ray flare) . 

For a clearer representation, the RXTE/ASM 
and Swift/BAT night-by-night data points in Fig- 
ures[T]to[3]were combined until one of the following 
conditions was met: (i) the combined data point 
had a statistical significance of more than 3 stan- 
dard deviations (cr), or (ii) 15 bins of the original 
light curve were combined. The RXTE/ASM and 
Swift/BAT data in Figure [4] are shown in dwell- 
by-dwell bins. The TeV data (Whipple and VER- 
ITAS) in all figures are shown in a run-by-run bin- 
ning of 10 — 20 min duration. All other data are 
shown in a binning corresponding to the point- 
ings/exposures of the individual experiments. 

Although the X-ray fluxes show a long-term 
structure with phases of higher activity followed 
by phases of lower activity, flux variations by a 
factor of two can be observed on time scales down 
to a few days. In the optical band, flux variations 
are observed on longer time scales of the order of 
a weeks to months. The structure of the optical 
light curve (i.e. Figure [Hand possible connections 
to the X-ray/TeV band are discussed in Sec. 13.41 
No significant flux variations can be observed at 
radio energies. The flux correlations between the 
different energy bands are discussed in Section l3~4l 
Two phases of enhanced X-ray and/or TeV activ- 
ity can be identified in the light curve: 

• Phase 1: The first active phase (phase 1) 
occurred during the summer 2006 and lasted 
for at least a few months. The active state 
can be identified in soft and hard X-rays 
(Figure[lj). During this time period no VER- 
ITAS data were taken and only a few nights 
,,„ , , nl 77 , . . , are covered by Whipple data. 

ii Mrk421 on the MOJAVE (monitoring of jets m active J ^ 

galactic nuclei with VLBA experiments) project page: m . . , , , 

|http://www.physics.purdue.edu /astro/MOJAVE/sourcepages/ nrri+^ shMt aSe 2 The SeCOnd actlve P^f e (P haSe 2 ) 

(TcTTurred in April/May 



Aller et al. I (|1985t ) . Some daily observations were 
averaged to improve the signal-to-noise ratio. 

Unfortunately, the source is rather weak which 
may mask some variability. Nevertheless, the UM- 
RAO measurements over many decades have iden- 
tified continuous fluctuations in amplitude but 
not a single outburst-like flare. Small structural 
changes in the radio jet are apparent in the MO- 
JAVE 15 GHz VLBI images for the soured^- 

3. Results 

In the whole VERITAS data set an excess of 
29974 7-ray events was detected from the direc- 
tion of Mrk421 after application of event selection 
cuts (31523 ON events, 7746 OFF events, normal- 
ization a = 0.2), corresponding to a statistical sig- 
nificance of 277 standard deviations. An overview 
of the light curves at radio to TeV energies is given 
in Section l3~Tl Subsequently, we discuss the time 
and spectral variability of the TeV 7-ray data (Sec- 
tion l3.2[) as well as the time and spectral variability 
of the X-ray fluxes (Section 13. 3p on different time 
scales. Finally, we scrutinize how those relate to 



2008 and was 



9 



recorded with excellent coverage in the X- 
ray and TeV bands (Figure [2]). However, 
a zoom-in ol this second flaring phase (Fig- 
ures El and 0| shows that the strongest 
TeV emission (phase 2b) is not coincident 
with the strongest soft/hard X-ray activity 
(phase 2a, peaking roughly one month before 
phase 2b) . The lack of increased X-ray emis- 
sion during the peak T eV flaring might in- 



dicat e an orphan flare (jKrawczvnski et al 
20041 ) . However, the characteristic time 



3.2. 



scales of flux changes in the TeV band can 
be less than an hour (the major flare is 
fully contained within a time interval of 
5 h) , so that a detailed comparison has to 
be restricted to closely-simultaneous data, 
see Section EH The TeV flare is followed 
by a somewhat enhanced X-ray flux: The 
Swift/BAT, Swift /XRT and RXTE/PCA 
data indicate a doubling in flux level between 
the night of the flare and the following night, 
declining back to the previous level within 
a few days, which is nicely sampled by the 
Suzaku/XIS (Figured]). The corresponding 
structure of the X-ray light curve, however, 
does not substantially differ from low-state 
variations, so that a physical connection to 
the TeV activity cannot be claimed. 

Temporal and Spectral Variability in 
the TeV 7-Ray Band 



Flux variability. Flux variations on time scales 
of 1-2 days are found in the TeV band. Except for 
two nights measured during a strong flare in May 
2008 (phase 2b, Figure EJ, no significant TeV flux 
variations are observed within individual nights. 
It should however be mentioned, that the obser- 
vation time during individual nights often did not 
exceed 10 — 20 mirl 12 !. This prevents placing strong 
constraints on the >0.5— 1 h time scale variabil- 
ity, given the low/medium states of Mrk421 dur- 
ing most of the measurements. Nevertheless, the 
strong outburst measured in May 2008 (Figure |4]) 
clearly shows variability on sub-hour time scales 
even though the flare was recorded at zenith an- 
gles down to 60 deg (see inserts in Figure SJ . At 



high zenith angles z, the effective areas vary dras- 
tically with a small change in z and the sensi- 
tivity suffers from the strongly increased energy 
threshold (leading to a huge loss in event statis- 
tics). Since the data points in the light curve are 
given above 300 GeV, the calculated fluxes (for the 
zenith angle range of 40 deg —60 deg) are derived 
(extrapolated) based on the measured spectral 
shapes (dN/dE = I -(E/1 TeV)- r -exp(-E/E cut )). 
Therefore, an increased systematic error on the 
integral flux of —30% for z > 45 deg and ~40% 
for z > 55 deg is assumed. Given these facts, we 
do not determine a quantitative value for the flux 
doubling times observed during this flare. 

Spectral variability. To investigate a possible 
change of the spectral shape as a function of flux 
state, the data are divided into subsets according 
to different flux levels. The separation into flux 
intervals is chosen such that reasonable statistics 
are guaranteed for each subset. An energy spec- 
trum is derived (see Section |2~T1) for each subset 
and is subsequently fit by a power-law function 
with exponential cutoff dN/dE = I -(E/1 TeV)~ r - 
exp(— E/E cu t). A fit of a simple power-law func- 
tion can be excluded with high confidence for most 
of the spectra (see as an example the dotted line in 
Fig. El resulting in a x 2 /dof = 264.3/15). The re- 
sults of the fits are summarized in Table [2] and the 
energy spectra are shown in Figure [5j No correla- 
tion between the cutoff energy E cut and the flux 
normalization Iq can be claimed. The spectra are 
also fitted with the same function by fixing the cut- 
off energy to S cut = 4 TeV (Table EJ. In this case, 
a hardening of the spectrum with increasing flux 
level can be seen, see Figure El A linear correla- 
tion between the flux and the index T is disfavored 
(x 2 /dof = 44.6/5, p = 1.7 • 10~ 08 ) as compared to 
a quadratic relationship T(Iq) = a + b ■ Iq + c ■ Iq 
with a = 2.70 ± 0.04, b = (-5.3 ± 0.6) • 10 9 and 
c = (1.08±0.20)-10 19 (x 2 /dof = 13.2/4, p= 0.01). 
This finding indicates that the spectral hardening 
with flux level flattens at very high flux values, 
as was seen already in th e case of PKS 2155-304 



([Aharonian et al. Il2009bl ). Given the sparse sam 



12 Regular run durations are 20min, but for monitoring pur- 
poses, some observations were conducted with just lOmin 
per night. 



pling during most of the nights, we were not able 
to further separate the data into rising and falling 
(with time) flux states which may have an effect 
on the r(io) function. However, the general flux 
versus T trend is in good agreement with earlier 
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Fig. 1. — Mrk421 light curves measured by different experiments (see Section [5]) in 2006-2008. All errors 
are statistical errors only. Shown are the radio band (upper panel), the optical band (second panel), the 
X-ray band (third & fourth panel) and the TeV band (bottom). Some of the RXTE/ ASM and Swift/BAT 
flux points have been re-binned for better visibility. Two phases of activity can be identified in the X-ray 
and TeV bands (phases 1 and 2). The vertical dotted line indicates the maximum of the hard X-ray flux 
(Swift/BAT, phase 1). During phase 2 a very good X-ray/TeV coverage was achieved, where the X-ray 
observations were partly triggered by VERITAS. 
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Fig. 2. — 2008 Mrk421 light curves. For details see caption and legends in FigureQ] The vertical lines enclose 
phases of enhanced X-ray (phase 2a) and enhanced TeV (phase 2b) activity; see Figure [3] for a zoomed view 
of these phases. An exponential rise/fall function was fitted to the optical R band Tuorla/KVA data ('Opt 1' 
& 'Opt2', see discussion in Sec. 13.41 and fit results in Tab. [3]). The dashed line represents the optical rise/fall 
function folded with the ASM X-ray (dwell- by-dwell) light curve (see text for disussion). The vertical arrows 
pointing downwards in the bottom panel indicate the dates for which we compiled quasi-simultaneous SEDs 
and modeled them with a SSC code (see Figures ITT1 fc [T^)) . 
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Fig. 3. — Mrk421 light curves measured during the X-ray/TeV flaring state in April/May 2008 (phase 2). 
For details see caption and legends in Figure Q] The vertical dotted line highlights the brightest X-ray flare 
observed with the RXTE/PCA. The vertical arrows pointing downwards in the bottom panel indicate the 
dates for which we compiled quasi-simultaneous SEDs and modeled them with a SSC code (Figures [TT] & 
CE2). 
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Fig. 4. — Mrk421 light curves measured during the TeV flaring state in April/May 2008 (phase 2b). For 
details see caption and legends in Figure [TJ Note, that the highest TeV flux around MJD 54590 was not 
accompanied by an increase in X-ray activity (compare with Figure fT0|) . The continuous Suzaku/XIS soft 
X-ray monitoring (starting two days after the strong TeV flare) is shown, as well. The inserts in the lowest 
panel show the VERITAS TeV light curves in a 3 min binning, covering the flaring time span indicated by 
the colored, vertical boxes. The flares were recorded under mostly large zenith angles (indicated by the 
arrows) with a correspondingly higher systematic error on the flux, see text for more details. 




Energy [TeV] 



Fig. 5. — Time-averaged VERITAS energy spectra for different flux levels. Only data taken in 2008 are used. 
A power-law with exponential cutoff is fit to each spectrum (curves), the fit parameters and statistical errors 
are summarized in Table [2] The dotted line shows a power-law fit to the low-state spectrum to illustrate the 
incompatibility of this model and the data. The high-flux spectra were taken at zenith angles z > 35 deg 
with a correspondingly higher energy threshold. 



results obtained with th e Whipple 10 m telescope 



( Krennrich et al. 1 120021 ) which are also shown in 



Figure [5] 

Mrk421 is detected by VERITAS (on average) 
with a statistical significance well above 10 stan- 
dard deviations per run, even for runs with a du- 
ration of only 10 min. This enables a run-by-run 
derivation of the energy spectrum on time inter- 
vals of 10 min or less. Any energy spectrum de- 
rived from an individual data run which meets 
the following requirements is fit by a power-law 
with exponential cutoff (the cutoff energy again 
being fixed to E cut = 4 TeV): (i) A differential 
flux point is only considered in a fit if the sta- 
tistical significance of the excess is above 2 stan- 
dard deviations; (ii) An energy spectrum is only 
fit if at least four differential flux points fulfill the 



first criterion. The results of the energy spectra 
derived for the individual runs are shown in Fig- 
ure [6] (gray points) and confirm the trend which 
has been found already in the data sets divided 
according to the different flux levels. 

For both cases, the correlation coefficient was 
calculated and the corresponding non-directional 
chance probability for the null hypothesis p (non- 
correlation) was calculated using the Student t- 
distributiorF^I. The correlation coefficients are r = 
—0.86 for the data sets separated by flux levef^l. 
and r = —0.55 (p < 10~ 4 ) for the distribution 
based on the run-by-run data sets. 



13 Note: The correlation factor does not account for the sta- 
tistical errors on the individual data points. 

14 Since this sample consists of only 7 data pairs no chance 
probability p was calculated for r. 
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Table 2: Parameters and statistical errors of fits dN/dE = Io • (E/lTeV)~ r ■ exp(— E/E cut ) to the energy 
spectra shown in Figure [5] observation time Xh vo , statistical significance of the excess in the data set, flux 
normalization Iq, photon index T, cut-off energy E cu t and the quality of the fit. The second row for each 
flux state corresponds to a fit in which the cut-off energy was fixed to a value of E cu t = 4 TcV. 





Tlive[h] 


sign, [a] 


10 11 

Io [ C m- 3 s- 1 TcV- 1 ] 




r 


£ C ut [TeV] 


X 


2 /dof 


very-low-state 


7.83 


77.62 


4.78 ±0.73 




2.29 ±0.11 


1.59 ±0.34 


17.8/11 


(1.62) 








3.056 ±0.086 


2 


608 ± 0.033 


4±0 


29.2/12 


(2.44) 


low-state 


16.3 


181.3 


7.60 ±0.33 


2 


285 ±0.035 


2.95 ±0.29 


20.7/14 


(1.48) 








6.769 ±0.084 


2 


375 ±0.015 


4±0 


29.3/15 


(1.95) 


mid-state 


7.79 


160.6 


10.26 ±0.44 


2 


278 ±0.037 


4.36 ±0.58 


14.9/15 


(0.99) 








10.54 ±0.15 


2 


256 ±0.017 


4±0 


15.4/16 


(0.96) 


high-state A 


1.38 


73.0 


19.08 ±2.63 




2.01 ±0.12 


1.91 ±0.41 


15.6/12 


(1.30) 








13.53 ±0.38 


2 


295 ±0.037 


4±0 


24.0/13 


(1.84) 


high-state B 


0.63 


47.3 


22.23 ±2.42 




1.88 ±0.24 


3.06 ± 1.07 


13.4/6 


(2.24) 








20.77 ± 1.14 


2 


053 ± 0.080 


4±0 


143.0/7 


(2.00) 


high-state C 


0.63 


53.76 


21.74 ± 1.97 




2.40 ±0.26 


9.6 ±9.0 


0.9/6 


(0.15) 








23.96 ± 1.62 




2.05 ±0.09 


4±0 


2.7/7 


(0.39) 


very-high-state 


0.63 


63.5 


35.77 ±3.08 




1.87 ±0.17 


2.74 ±0.60 


6.2/9 


(0.69) 








32.0 ± 1.2 


2 


111 ±0.057 


4±0 


8.7/10 


(0.87) 



3.3. Temporal and Spectral Variability in 
the X-ray Band 

Energy spectra and fluxes. The RXTE/PCA 
spectra were fitted in an energy range of 3— 20 keV, 
while the Swift/XRT spectra were fitted between 
0.4 — lOkeV. Two models were tested to fit 
the data: a power law $>(E) — kE~ T and a 
log-parabolic model p(ff) = fc.fiM 0+61o s(-E/Eo)) 
(|Massaro et al. I I2004L \20m . The log-parabolic 



function uses an energy dependent photon index 
T(E) = a + blog(E / Eo); the parameter b defines 
the curvature in the logarithmic parabola, and a 
is the spectral index at Eq. Both models account 
for absorption assuming a fi xed galactic column 



densi ty of 1.61 x 10 cm (jLockman fc Savage 



1995). The mean reduced x values from log- 
parabolic fits to Swift/XRT and RXTE/PCA data 
are 1.20 (96 dof) and 0.87 (161 dof), respectively. 
These are significantly lower than the respective 
mean reduced \ 2 values from power law fits of 
1.80 and 1.35. Hence, the integral flux from 
2—10 keV was calculated for each observation from 
log parabolic fits. 



Flux variability and light curves. The X-ray 
light curves of the 2 — lOkeV fluxes are shown in 
Figures Q] to [4] Significant day time-scale flaring 
is seen in many months in both 2006 and 2008, 
while for the observations in 2007 the Mrk421 
X-ray flux remained relatively low. Of particu- 
lar interest are the observations in March to May 
2008. RXTE/PCA recorded one of the highest 
ever X-ray fluxes for Mrk421 on March 30, 2008 
(phase 2a). However, during the very high TeV 
7-ray flux state measured by VERITAS in early 
May 2008 (phase 2b) the near-simultaneous X-ray 
flux was only found to be at a moderate level. 

For all but one case the RXTE observation ex- 
posures were shorter than 1.5 hours (with an av- 
erage exposure of 25 min). The Swift observations 
had a far larger range in exposures, from 15 min- 
utes in many observations, to over 40 hours in 
each of seven observations; the average exposure 
per pointing was 3.3 hrs. Significant variability 
in both RXTE/PCA and Swift/XRT observations 
was found. The corresponding Swift observations 
spanned up to 3 days, with clear flaring in the rates 
on hour time-scales. Some of the RXTE/PCA 
light curves show a steady change in the count 
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Fig. 6. — Photon index T versus flux normaliza- 
tion Jo (at 1 TeV) obtained from a power-law fit 
with exponential cutoff (fixed to E'cut = 4 TeV, 
statistical errors only). Each gray data point cor- 
responds to a data run with a duration of 10 to 
20min. The black points correspond to the val- 
ues obtained from the fits to the spectra grouped 
according to their flux states (defined in Figure [5] 
and Table [5]). Correlation coefficients r and the 
corresponding chance probabilities for the null hy- 
pothesis (no correlation) are given in the legend. 
Results obtained from earlier Whipple 10m obser- 



2002) are shown for ref- 



vations (jKrennrich et al. 
erence. 



rate over 20 minute periods. However, more de- 
tailed studies of the sub-day flux variations in the 
X-ray data are beyond the scope of this paper. 

Spectral variability. Drawing from this large 
set of X-ray observations, the correlation of X- 
ray flux to spectral shape is investigated here. 
This will allow identification of trends which could 
give important input to the modeling. The left 
panel of Figure [7] shows the correlation plot be- 
tween the log-parabola parameters a (index) vs. 
the 2 — lOkeV flux. The corresponding corre- 
lation plot between the curvature parameter b 
and the flux is shown in the right panel of Fig- 
ure [7] Although the index a is not independent of 
the curvature parameter b, the conclusion can be 
reached that an increased flux is accompanied by 
a hardening of the spectrum (parameter a) since 
there is only a moderate correlation between the 



flux and b. This finding is compatible w ith ear- 
lier findings (jFidelis fc Iakubovskvil 120081 ). The 
(anti)correlation between a and the flux is sig- 
nificant: the non-directional chance probability 
of the derived correlation coefficient for the null- 
hypothesis is p < 10~ 4 for the RXTE/PCA data, 
as well as for the Swift/XRT data. 

3.4. Flux Correlation Analysis 

This section describes the search for flux corre- 
lations between the light curves measured in differ- 
ent energy bands and is based on the light curves 
shown in Figures [TJ to |U 

Radio/optical/ TeV flux correlations. As 

can be seen in Figure Q] there is no significant vari- 
ation in the radio flux on day, week, or month 
time scales. Thus, no significant correlation with 
fluxes measured in the other wave bands is found. 
This could be explained if (i) the dominant portion 
of the radio emission does not originate from the 
inner jet region that likely produces the flux vari- 
ability in the other energy bands, or (ii) structures 
in the light curve are smeared out due to slower 
cooling of the radio emitting electrons compared 
to the electrons responsible for the shorter wave- 
length emission (see Sec. [4j. 

The optical fluxes, on the other hand, show 
clear variability on time scales of weeks (Figure [2]), 
where the fluxes from the different bands (R, B, 
and V) are clearly correlated. This is not surpris- 
ing since the frequency filter bands are not largely 
separated in terms of photon energy and a com- 
mon origin of the radiation in the three bands can 
be assumed. Figure [8] shows the correlation plot 



Table 3: Parameters of an exponential rise/fall 
function (Sec. O $ op t(*) = a + b/{e-^- T ^' T ' + 
e (t-T D )/ Tf j fitted t0 the pti C al flares 'Opt 1' and 

'Opt 2' (Fig. EJ. 



Parameter 



Optl 



Opt 2 



a [Jy] 0.021 ± 0.001 0.023 ± 0.001 

b [Jy] 0.023 ± 0.002 0.024 ±0.003 

r r [d] 7.16 ± 1.78 9.52 ±1.35 

T f [d] 5.46 ± 0.80 4.46 ± 0.77 

Tp [MJD] 54476.8 ± 1.4 54501.4 ±0.7 
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Fig. 7. — Correlations between the spectral shape and the flux level for the X-ray data sets taken with 
RXTE/PCA and Swift/XRT. Correlation coefficients r and the corresponding chance probability for the 
null hypothesis (no correlation) are given in the legends. Left: The log-parabola parameter a vs. the 
2 — lOkeV photon flux. Right: The log-parabola curvature parameter b vs. the 2 — lOkeV photon flux. 



between the VERITAS TeV 7-ray fluxes and the 
optical fluxes in the R, B, and V bands. A maxi- 
mum time gap between the TeV and optical mea- 
surements of At < 0.5 days was allowed for the 
individual data pairs. The corresponding correla- 
tion factors are compatible with chance expecta- 
tions. A similar correlation study between X-ray 
and optical fluxes did not result in any significant 
correlation coefficient, either. Given the different 
time scales of the flux variations (sub-day level in 
the X-ray and TeV band, and weeks in the opti- 
cal band), the lack of a direct flux correlation is 
expected. 

Except for the Swift /UVOT data the contribu- 
tion of the host galaxy is not subtracted. Since 
it has to be constant in time, the light curves can 
be used to set an upper limit on the base-line, 
which in turn would be an upper limit on the host 
galaxy contribution in the measured data. This 
contribution is estimated to be 0.009 Jy (R band, 
see dotted line in Fig. [Q), 0.009 Jy (V band), and 
0.006 Jy (B band). These estimates are somewhat 
lower as compared the modeling of iNilsson et al. 
(|l999h who estimate 0.014 Jy for the R band. 

Although neither instantaneous nor delayed 



correlation between the optical and the X-ray/TeV 
bands was found it is important to study the struc- 
ture of the optical light curves. Two well defined 
flares^! 'Optl' and 'Opt 2' (Figure [2} were fitted 
with an exponential rise/fall function & pt(t) = 
a + 6/(exp(-^p) + exp(^p)). The fit param- 
eters are summarized in Tab. [31 showing that the 
optical flux changes on time scales of less than 
ten days with the indication of slightly shorter fall 
times Tf as compared to the rise times r r . 

These time scales are on the same order as the 
estimated value of ~15d, based on the character- 
istic X-ray/TeV variability time scale of < 0.5 d 
(see the discussion in Sec. HJ). One can inter- 
pret <f?opt(i) as a 'characteristic' optical response 
to a single X-ray/TeV flare. With this assump- 
tion a hypothetical prediction can be made for 
the optical light curve by folding $ op t(i) (using 
r r = 7d and Tf = 5d, and a delay AT = 7d) 
with the X-ray and/or TeV 7-ray light curves. 
The prediction based on the frequently sampled 
RXTE/ASM X-ray light curve (dwell-by-dwell, 



A third flare occuring around MJD 54149 was also fitted, 
resulting in comparable structural properties as listed in 
Tab. [3] 
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Fig. 8. — Correlation between TeV and optical 
fluxes. The maximum allowed time difference be- 
tween the data points is At < 0.5 days. 

MJD 54450 to 54630, linear interpolation between 
the flux points) is shown together with the op- 
tical light curve in Figure [2] (arbitrary units) - 
no agreement is found. A caveat should be men- 
tioned: the RXTE/ASM flux measurements are 
not very accurate; however, the general trends 
(high- vs. low-state) should allow the compari- 
son of the two different wave bands. $ pt(£) was 
also folded with the TeV 7-ray light curve (MJD 
54500 to 54630, linear interpolation between the 
flux points). Again, no correlation can be found. 
There is a caveat here, as well: the non-continuous 
sampling of the TeV 7-ray light curve and the 
short duty-cycle of flares results in a considerable 
chance that one or more strong 7-ray flares have 
been missed which would change the shape of the 
folded optical light curve. Therefore, the latter 
results are not shown in Figure [2] Furthermore, 
our above comparison ignores the possibility that 
particles with different energies are injected at the 
beginning of the flare. 

Hard/soft X-ray flux correlations. The cor- 
relation between hard and soft X-ray fluxes is 
studied based on the dwell-by-dwell flux points 
measured with Swift/BAT (15 - 150 keV) and 
RXTE/ASM (2 - 10 keV) . The correlation plot in- 



Fig. 9. — Hard/soft X-ray flux correlation be- 
tween the 15 - 150 keV Swift/BAT flux and the 
2 - lOkeV RXTE/ASM flux. All dwell-by-dwell 
flux points with a statistical significance above 2 a 
measured in the first half of 2008 (MJD 54457- 
54640) are shown for two cases of a maximum time 
lag between the observed data points. The legend 
gives the corresponding correlation coefficient and 
the corresponding chance probability for the non- 
correlation. 



eluding all data points taken in spring 2008 (MJD 
54457-54640, contemporaneous with the VERF 
TAS coverage) with a statistical significance of 
more than 2 standard deviations is shown in Fig- 
ure [9j The distribution was generated for two 
different requirements regarding the maximum al- 
lowed time gap between the center time of the in- 
dividual pointings (At < lh and At < lOmin). In 
both cases, indications for only a weak correlation 
of r ~ 0.1 are found with moderate significance. 
However, while the RXTE/ASM data are testing 
the falling edge of the synchrotron peak in the 
SED, the Swift/BAT data likely fall into the tran- 
sition zone between synchrotron and high-energy 
peaks, compare with Figures [TT1 and IT2"1 

X-ray/TeV flux correlations. The rich data 
sets allow for a detailed study of the X-ray versus 
TeV flux correlations. The TeV runs have a dura- 
tion of 10 — 20 min. With a few exceptions (during 
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Fig. 10. — Correlation between TeV fluxes (VER- 
ITAS and Whipple) and X-ray fluxes (RXTE: 
PCA & ASM, Swift: XRT & BAT). All X-ray 
fluxes are given between 2 — 10 keV, except for 
the Swift/BAT fluxes which are given between 
15 — 150 keV. Only data measured in the first 
half of 2008 (MJD 54457-54640) are shown. The 
RXTE/ ASM and Swift/BAT fluxes are dwell-by- 
dwell with the requirement that the statistical sig- 
nificance is above 2 standard deviations. The max- 
imum time difference between the data points is 
At < 15min. The correlation coefficients includ- 
ing the chance probability for the null hypothesis 
are shown in the legend. The flux pairs measured 
during the first (2b(l)) and second night (2b(2)) 
of phase 2b are marked. 



the highest flare) no indications for TeV flux vari- 
ations were found within any of the runs. A max- 
imum time lag between the correlated X-ray/TeV 
data points of At < 15min was allowed (close 
to the average TeV run duration). Since some of 
the dedicated, high-quality data from Swift/XRT 
and RXTE/PCA were considerably longer in ex- 
posure, an additional cut on the total length of 
the X-ray pointing was applied (based on the av- 
erage exposure per pointing): AT < 30min in the 
case of RXTE/PCA and AT < 2h in the case 
of Swift/XRT. Flux variations within individual 
X-ray pointings are not investigated in the frame- 
work of this paper; therefore, the data pairs can- 



not be considered as exactly simultaneous. The 
RXTE/ ASM and Swift/BAT dwell-by-dwell data 
points have exposures of the order of 1 — 10 min, 
so that no cut on AT was applied in those cases. 
Individual RXTE/ ASM and Swift/BAT flux mea- 
surements were only considered in case of a signif- 
icance level of > 2 a. The results discussed below 
do not strongly depend on the exact choice of At 
and AT. 

Figure [TUJ shows the flux correlations for the 
data taken in spring 2008 (MJD 54457-54640). 
The correlation coefficients and the corresponding 
non-directional chance probability for the null hy- 
pothesis (no correlation) were calculated and are 
shown in the figure legend. The X-ray and TeV 
fluxes seem to be correlated (r > 0.5 for most data 
subsets) with chance probabilities of the order of 
a few percent or below. The VERITAS/BAT data 
set is the only one which is not correlated in a sig- 
nificant manner. One of the high TeV flux points 
measured during the strong TeV 7-ray flare (Fig- 
ure I3J was accompanied by a Swift/BAT dwell- 
by-dwell pointing (At < 3.4 min), which however 
did not indicate an increased activity in the hard 
X-ray band. The Swift/BAT point has a dura- 
tion of 8 min and is therefore fully contained in 
the time intervall of the corresponding VERITAS 
data run of 20 min duration. This flux pair is lo- 
cated in the lower right corner of Figure [10] (la- 
beled as 2b(l)) and may be seen as the indication 
of an orphan TeV 7-ray flare. However, since this 
indication is based on only one X-ray/TeV flux 
pair measured during a high TeV 7-ray flux state, 
there is not enough evidence for a strong claim. 
If the corresponding data pair were removed, the 
VERITAS/BAT flux correlation would increase to 
r = 0.62 with a chance probability of p — 0.019. 
Other flux pairs from the two TeV flare nights 
(phase 2b) are indicated in Figure [TU] as well. 

Discrete correlation functions. The light 
curves were also analyzed using the dis- 
crete correlation funct ion (DCF) technique 
(|Edelson fc Krolik I Il988t ) in order to search for 
correlated/delayed emission between different en- 
ergy bands, allowing for significant time lags (i.e. 
one light curve having the delayed shape of an- 
other one). The radio, optical and X-ray light 
curves were tested against the TeV 7-ray light 
curve. Except for the zero lag X-ray /TeV cor- 
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relation (see above) no significant time lag and/or 
correlation at zero time lag was found for the ra- 
dio and optical band as compared to the TeV 7-ray 
emission. 

3.5. Spectral Energy Distributions and 
Modeling 

The spectral energy distributions (SEDs) show- 
ing the VERITAS and multi-wavelength (MWL) 
data were generated for individual nights and plot- 
ted in Figures Qj] and Q21 Based on the differ- 
ent variability time scales, data from the different 
wavebands are plotted in the quasi-simultaneous 
SEDs if the time lag At between the MWL data 
and the TeV data is Aix-ray < 0.15 d (X-ray), 
At opt < 1.5 d (optical), and At la< a < 3d (radio). 
Shorter time lags are present (down to real simul- 
taneity) and the exact times are given in the figure 
legends of the SEDs. The time spans (duration) 
of the X-ray observations are given in the legends, 
as well. Possible spectral variations within indi- 
vidual X-ray observations are beyond the scope 
of this paper and are ignored in the modeling of 
the SEDs. The X-ray spectra (Swift/XRT and 
RXTE/PCA) were fit with a log-parabola model 
= feJ5-(°+ 61o s(-B)) ; an d the corresponding 
log-parabola bow ties including the statistical er- 
rors of the fits are shown if the fit resulted in a 
X 2 /dof < 2. All errors are statistical only. 

As a reference, Figures [TT] and Q21 also show 
the non-simultaneous MeV/GeV energy spec- 



trum measured by the Fermi/LAT (Atw ood et al. 



2009 ) in August 2008, to July 2009 (|Abdo et al 

2010) . During this period, the LAT detects 



Mrk421 with a statistical significance of 88 a and 
the energy spectrum is well fit by a power law 
in the whole energy range (curvature index given 
in the Fermi catalog of 0.72). Interestingly, the 
variability index of 43.9 does not indicate very 
strong variability during the above period, so 
that the Fermi/LAT energy spectrum shown in 
the SEDs may give a realistic indication for a 
low/medium flux state of Mrk421. However, given 
the non-simultaneity of the Fermi observations, 
the MeV/GeV spectrum was not included in the 
modeling. 

The S EDs were fit with a one-zone SSC model 
following iKrawczvnski et al. 



(|2004h . 



In this 

model a spherical emission region of radius R is 
filled with a relativistic electron population, trav- 



eling down the jet with a bulk Lorentz factor T 
((3 = v/c). The magnetic field B in the emission 
region is randomly oriented. The emitted radia- 
tion is Dopplcr-shiftcd by 8 = [T(l — /3cos0)] _1 , 
9 being the angle between the jet axis and the 
line of sight of the observer. The electron distri- 
bution is normalized by a factor u e (in units of 
ergs/cm 3 ) and is described in the jet frame by a 
broken power law with E m i n , -Ebrk, and E max and 
the two corresponding spectral indexes a% — —2.2 
and ct2 — —3.2. Mrk421 is located at a redshift of 
0.031 so that the energy-dependent pair absorp- 
tion on photons of the extragalactic background 
light (EBL) cannot be neglected. EBL absorp- 
tion is taken into account in the model following 
Franceschini et al. (2008). Since most of the op- 
tical data were not corrected for the contribution 
of the host galaxy, the optical data points have 
to be seen as an upper limit on the jet emission 
region. The estimate on the host galaxy contri- 
bution (Sec. 13. 4p has been added as a systematic 
error. The SSC model was adjusted to the SEDs 
with the following procedure: 

• In a first step, the three SEDs with the 
smallest time lags between the X-ray and 
the TeV spectra (MJD 54559.25, 54562.25, 
and 54509.37) were used to derive a full set 
of model parameters. The Doppler factor 
was set to 8 = 40 and the magnetic field 
was set to B = 0.2 G, in agreement (order 
of magnitude) with a synchrotron cooling 
time of T sync (B) x 8 rj Tfl are = lh for the 
electrons emitting the synchrotron emission 
close to the maximum of the spectral energy 
distribution. Then, the radius R, the energy 
density in electrons u e [ergs/cm 3 ], as well as 
-Ebrk and E max were varied until the model 
would describe the SED (optical to TeV). 

• In a second step, the remaining SEDs were 
fit by only allowing variations of u c and 
-Emax (and in case of bad fits also Bbrk) as 
compared to the model parameters derived 
above, reflecting a change in the injected 
electron population. 

The SSC model fits are shown in Figure HT1 and 
1121 and the model parameters for the individual 
SEDs are summarized in Table |4j The models 
generally under-predict the radio emission which 
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is synchrotron self-absorbed in the model at the 
low-energy tail of the SED. This discrepancy could 
be explained by additional radio emission from re- 
gions in the jet not emitting the TeV radiation. 
Also, the synchrotron self-absorbed radio blobs 
could expand and le ad to a delayed radio e mission 



from a larger region (jAcciari et al. Il2009al ) , which 



however is not taken into account in the model. 

The break energy E^k of the electron spec- 
trum is very near the maximum energy of the 
electron spectrum E max . The missing plateau 
(Ebrk ~ Bmax) is directly seen in the Swift/XRT 
data which show a direct turn-over of the syn- 
chrotron peak. This is in agree ment with ear 



lier fi ndings in the case of Mrk421 (Fos sati et al. 



2008) and is in contrast to the SEDs measured 
in the case of Mrk501 for which the synchrotron 
emission peaks at hig her energies and shows in - 



2000). 



dications of a plateau ([Krawczvnski et al. 
If E m - m <C -Ebrk, the missing plateau implies that 
cooling did not have time to kick in. However, 
if -E m m = ^brk and cooling is very efficient, elec- 
trons cool below -Ebrk and one gets a new 'effective' 
E m i n <C Ebrk, and the 'true' E m - ln becomes Ebrk- 

All SEDs are reasonably described by using the 
same model parameters except for the electron 
normalization u c and the break/maximum elec- 
tron energies Eb r k and E max which were adjusted 
for each individual SED (Table 0). The only ex- 
ceptions are the X-ray flare (Figure ITT]) and the 
two TeV 7-ray flare SEDs (Figure [JJ) : 

• MJD 54555.38: The RXTE/PCA bowtie 
shown in this SED (Figure [TTj) does not 
qualify for the previously defined X-ray/TeV 
time lag criterion. However, the X-ray flux is 
the highest one measured during the whole 
MWL campaign (see dotted line in Figure |3]) 
so that it is shown for reference (but not in- 
cluded in the fit). The closest Whipple TeV 
flux point is ^2 h away and also does not 
show signs for an increased TeV flux. This 
may indicate an orphan X-ray flare, but the 
lack of simultaneous TeV data does not al- 
low a strong conclusion. 

• MJD 54588.3: The SSC model slightly over- 
predicts the optical emission. However, the 
optical data are not simultaneous and the 
TeV fluxes can change within 20 min: the 
black data points in Figure [T5] represent the 



first 20 min of the flare, whereas the open 
gray points represent the second 20 min of 
the flare in that night (compare with the 
lower left insert in Figure 01 . Unfortunately, 
no X-ray data (Swift/XRT or RXTE/PCA) 
were taken during this night. 

• MJD 54589.21: The black data points in 
Figure [T5] represent the low plateau at the 
beginning of the flare - compare with the 
lower right insert in Figure |4] - and are 
the ones which were used for the fit of the 
SED. The RXTE/PCA measurement par- 
tially overlapped this time, but is over- 
predicted by the model. Given the gener- 
ally good fits of the model, this is interest- 
ing since the VERITAS/BAT flux pairs of 
the same night (see 2b(2) in Figure fTUj) as 
well as the previous night (2b(l)) also in- 
dicate a possible orphan TeV 7-ray flare. 
The higher flux states during this flare night 
(open gray points in Figure IT2|) were not ac- 
companied by simultaneous X-ray measure- 
ments and were therefore not modeled. 

The fact that the X-ray (and optical) data are 
not perfectly described during the flare nights (as 
well as for a few other nights) may be explained 
by the fact that only u c , Ebrk and E max were al- 
lowed to vary, after the other parameters had been 
fixed based on the three most complete and con- 
temporaneous SEDs, which all correspond to TeV 
low/medium states of Mrk421. 

4. Summary and Discussion 

Together with many MWL partners VERI- 
TAS conducted an intensive MWL campaign on 
Mrk421 in 2008. During low states of Mrk421, 
VERITAS is able to measure the source flux 
to an accuracy of 5% in 20 min time intervals. 
In the 47 hr data set presented in this paper, 
we did not find evidence for rapid flux variabil- 
ity on ti me scales of minu tes as reported for 
Mrk501 dAlbert et al. 1120071) and PKS 2155-304 



(Aharonian et al. 



2009bj T However, the only 
strong outburst (reaching a flux level of 10 Crab) 
which would have allowed testing these short time 
scales of flux variations was measured at high 
zenith angles with a strongly reduced detection 
sensitivity. 
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Fig. 11. — SEDs of Mrk421 for individual dates showing all available MWL data (statistical errors only) if 
the corresponding time lag criterion is met (see text for more detail). The legends indicate the time lags 
AT = Tmwl — TtcV between the individual MWL data and the VERITAS data. For the X-ray observations, 
the duration of the pointing A is given, as well. The number in parentheses refers to the phase (2a or 2b) 
during which the The SEDs were fit by a one-zone SSC model (dashed curves); the corresponding model 
parameters are summarized in Table 01 The non -simultaneous MeV /GeV spectrum measured by Fermi 
(August 2008 to July 2009) is shown for reference (|Abdo et al. Il2010h . 
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Fig. 12. — SEDs of Mrk421, continued from figure [TT] The SEDs of the two flare nights are also shown: 
MJD 54588.30 (no X-ray coverage in this night) and MJD 54589.21 with X-ray coverage during the onset of 
the flare (solid TeV points) but not during the high flare state (open TeV data points). 
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Table 4: Parameters of the SSC models as shown in 
Figure QT] and [T2] The following parameters are 
the same for all models: Doppler factor 6 = 40, 
magnetic field B = 2.0 x 10~ 5 T = 0.2 G, radius 
R = 2.5 x 10 15 cm, log(£: m i n ) = 3.0. The table 
shows the remaining parameters u c , log(£' max ) and 
log(£'b r k) for the different SEDs, as well as the 
particle energy density to magnetic field energy 
density ratio u e /uB = Sttu c /B 2 . 



date 




log(_E max ) 


log(£'brk) 


U e /UB 


[MJD] 


[ergs/cm 3 ] 








54475.4 


0.45 


10.8 


10.5 


283 


54478.4 


0.65 


10.9 


10.4 


408 


54508.3 


0.50 


11.0 


10.6 


314 


54509.3 
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11.1 


10.5 


377 


54536.4 
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251 
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ing time, one can estimate the expected relation 
between the observed time scales of flux varia- 
tions: The mean energy of emitted photons Sobs 
scales with the electrons' Lorentz factor squared, 
^obs oc 7 2 . The cooling time t coo i scales propor- 
tional to I/7. The energy of the X-ray photons 
is three orders of magnitude higher than the en- 
ergy of the optical photons. The Lorentz factors 
of the X-ray emitting electrons should thus be 1.5 
order of magnitudes (factor ^30) larger than the 
Lorentz factors of the optically emitting electrons. 
If the X-rays (and 7-rays, from inverse Compton 
scatterings of the electrons which emit X-rays as 
synchrotron emission ) show flux variability on a 
time scale of <, 0.5 d (|Fidelis fc Iakubovskvill20o3 
iGaidos et al. ~lll996l) . then the optical fluxes should 
vary on a time scale of < 15 d which is well com- 
patible with the observed time scales in the opti- 
cal and X-ray /TeV bands (Sec. 13. 4p . However, as 
discussed in Sec. 13.41 the optical light curve does 
not seem to be a delayed and stretched version of 
the X-ray light curve or the TeV 7-ray light curve, 
implying, that the dominant fraction of the ob- 
served optical emission does not originate f rom the 
X-ray /TeV emission region. Iterestingly, iRieger 



Two phases of enhanced X-ray (phase 1) and 
X-ray/TeV (phase 2) activity are found. Phase 2 
can be separated into a period of strong X-ray ac- 
tivity without a strong TeV 7-ray flaring, followed 
by a TeV 7-ray flare lasting for two days without 
any indication for contemporaneous strong X-ray 
activity. This may indicate an X-ray (phase 2a) 
and TeV 7-ray (phase 2b) orphan flare, but the 
data are too sparse for a definite claim. In the 
remaining data there is significant evidence for a 
correlation between the X-ray and TeV fluxes. 

No significant flux correlations between the TeV 
band and the optical/radio bands were found. As- 
suming that (i) the optical and X-ray/TeV pho- 
tons are emitted co-spatially, and that (ii) the 
flux variability time scale equals the radiative cool- 



(2004) discusses variable (even periodic) flux vari- 
ations on (periodic) time scales of P < 10 d which 
can be explained by geometrical arguments of in- 
ternal jet rotation. Our data, however, seems to 
indicate more complicated structures than period- 
icity. A similar estimate using the radiative cool- 
ing time for the energies in the radio band leads 
to time scales of 10 years or more, impossible to 
test with the given MWL data set. 

Clear indications for spectral hardening with 
increasing flux levels are found in the X-ray and 
TeV bands. In the TeV band, the spectral hard- 
ening seems to level out for the very high fluxes 
above ^5 Crab. A similar trend had already been 
found in earlier Whipple data of Mrk421, as well 
as in the str ong flare of PKS 2155-30 4 measured 



2009bD 



by H.E.S.S. (jAharonian et al 

The rich MWL data set presented in this 
paper allowed for the compilation of 18 quasi- 
simultaneous SEDs, well constrained by accurate 
Swift/XRT X-ray and VERITAS TeV spectra. 
The SEDs can be described by a one-zone SSC 
model with nearly one set of parameters; only u e , 
-E D rk and E max needed to be adjusted to describe 
the whole set of SEDs. However, for most of the 
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SEDs Mrk421 was found to be in low or moderate 
flux states, leading to similar SEDs on different 
days. Our SSC modeling indicates that the emis- 
sion process and the jet parameters are reasonably 
well constrained. 

So far, there are only four TeV blazars with a 
reasonable amount of simultaneous MWL data in 
order to claim a correlation between X-ray and 



TeV 7-ray fluxes: Mrk42 1 (IFossati et al. 



Blazeiowski et al. I 12005 



MrkSOI (Krawczynski et al. 



304 (Aharoni an et al. I 
2010h . 



(Acciari et al. 



Horan et al 
2002ft 



2009bJ) 



2008; 
2009), 



PKS 2155- 
and lES2344+5f2 
Such a correlation im- 
plies that the same high-energy particle pop- 
ulation (e.g. electrons) is responsible for the 
synchrotron emission at X-ray energies, as well 
as the high-energy IC emission at TeV ener- 
gies, as predicted in the framework of SSC mod- 
els. Investigating the exact shape of the X- 
ray/TeV correlation (linear, quadratic, etc.) will 
be one of the important goals for future stud- 
ies. Although this correlation is seen as a gen- 
eral trend, it does not nece ssarily hold true at 



the l evel of individual flares (|Krawczvnski et al 
2004). In our data, we find indications of an 



X-ray high-state not accompanied by TeV 7-ray 
flaring (phase 2a), as well as a TeV 7-ray flare 
without increased X-ray activity (phase 2b). Al- 
though the data is not exactly contemporane- 
ous - not allowing for a firm conclusion - such 
orphan flares in g eneral would require fine tuning 
of the SSC model (iKrawczvnski et al. I l2004h or al- 
ternative models, e.g. external-Compton models 
or models where the 7-ray emission is produced 
by hadrons, e.g. as proton - synchrotron emiss ion 
(|Miicke fc Protheroe I l200lt lAharonian 1 |2002|) or 



throu gh a proton induced cascade (jMannheim 
Il998 ). Further observations are needed to under- 



stand this particular aspect of the TeV flaring ac- 
tivity. 

reported comparable 
the Mrk42I X-ray 



Acciari et al. I (|2009bD 

week-scale trends between 



and optica l fluxes without a strong X-ray/TeV 



coupling. Aharonian et al. I ( 2009a ) reported 



clear indication of an optical/TeV flux correlation 
in case of PKS 2155-304. However, in the sec- 
ond PKS 2155-304 flare the optical/TeV correla- 
tion was not seen alth ough the data were agai n 
strictly simultaneous ( Aharonian et al.~l [2009b). 
Also, an optical/TeV correlation is not found in 



the large data sample presented in this paper or 
in earlier large data sets. Therefore, it does not 
seem to be a general property of TeV blazars. This 
may indicate that (part of) the optical emission is 
dominated by a region larger than the TeV emis- 
sion site and/or a different emission mechanism is 
at play. However, a certain level of optical syn- 
chrotron emission is un-avoidable given the syn- 
chrotron emission at X-ray energies. The ambigu- 
ous findings so far in terms of the optical/TeV cor- 
relation may further indicate that different emis- 
sion scenarios may play a role in different situ- 
ations, e.g. we do not always observe the same 
type of flares. No correlation between TeV and 
radio fluxes has been established at this point re- 
sulting in similar arguments as above concerning 
the emission regions/mechanisms. 

TeV flux variations are measured for dif- 
ferent TeV blazars with characteristic time 
scales down to 5 — 20min dAlbert et al. 



Aharonian et al. Il2009bl iGaidos et al. Ill996l) . In 



2007 



terestingly, the corresponding size of the emis- 
sion region reaches down to the order of 
the Schwarzschild radius o f the black hole of 
the corresponding AGN (jAlbert et al. I [2007; 



Aharonian et al. I l2009b| ). This can be seen as 



an indication that the TeV 7-ray emission from 
blazars comes from the base of the jet where the 
jet energy density is highest and the jet cross sec- 
tion is smallest. A similar finding was made in 
case of the radio galaxy M 87 for which a promis- 
ing approach to locate the site of the TeV emis- 
sion region was presented based on the combina- 
tion of TeV 7-ray observations wit h simultaneous 



high-r esolution radio observations ([Acciari et al. 
2009al) . 

We see two paths towards improving our 
understanding of the inner workings of AGN jets. 
The first path involves simultaneous TeV obser- 
vations with imaging telescopes with < 1 mas 
angula r resolution, e.g. VLBA ([Acciari et al. 
2009al ). and/or with polarimet r ic ob servations 



in the optical (jMarscher et al. I |2008) and the 
X-ray baneful. Furthermore, future observations 
in the MeV/GeV band with Fermi, together 
with observations in the GeV/TeV band with 
VERITAS, will facilitate the study of the spectral 
slope, time scales and correlation of the fluxes in 



3 see for example http://heasarc.gsfc.nasa.gov/docs/gerns/ 
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both energy bands which are important inputs 
for the theoretical modeling. The second path for 
achieving further progress concerns the theoretical 
modeling of the results. The analysis presented 
in this paper confirms that SSC models are 
successful in describing the broad-band spectral 
energy distributions of high-frequency peaked 
BL Lac objects. The two lessons which we 
infer from the modeling are (i) a large value of 
the relativistic Doppler factor 6 is required to 
explain the SEDs and the rapid flux varia bility 
(jGaidos et al. Ill996t [rTegelman et al~ll2008l) : (ii) 
in the case of SSC models, the electron energy 
density exceeds the magnetic field energy density 
both measured in the rest frame of the emitting 
volume (see Table SJ. In the case of EC models, 
one can find models with approximately equal 
electr on and magnetic fiel d ene rgy densities 
(e.g. Ghisellini fc Tavecchiol (<2009l )). or models 
in wh ich this ratio devia t es co nsiderably from 
unity iKrawczvnski et al. I (|2002h . Considering 



that protons may add to the particle energy 
density, the particle energy density may still 
be an equally comparable component in the 
plasma blobs producing HBL flares. Further 
progress can be achieved by combining them 
with those from theories describing the forma- 
tion and structure of jets. Recently, general 
relativistic 
codes 



(e.g. McKinncy (2006 


); Komissarov 


Krolik & Hawlevl ( 


2010); 


Spruitl (2010h 



have been used to validate aspects of analytic 
models of the magnetic form ation, accelera- 



1987 



tion, and collimation of jets (IWeber fc D avii 
19671: iBlandford fc Znaiekl Il977t iPhinnevI Il98 



Camenzind 19861: Lovelace. Wang fc Sulkanen 



Li. Chiueh fc B cgclman 




Vlahakis fc Konigll 12004 iKrawczvnskil 12007 
Giannios et al. 2009 discussed the model of 
"mini-jets in a jet" driven by magnetic re- 
connection. Mini-jets in jets may be able to 
reconcile the predictions of magnetic models 
of jet formation with the results from SSC (or 
EC) modeling of the blazar emission of jets: 
The reconnection mechanism conv erts mag- 



netic e nergy into particle energy (ISikora et al 
2005 : Giannios. Uzdenskv. fc Begelmanl .2009; 



Sikora et al.l 2009 ) and may thus be able to 
explain how particle dominated plasmas arc 
created in jets in magnetic field dominated 



jets. Furthermore, the creation of mini-jets 
with modest bulk Lorentz factors inside a jet 
with a modest bulk Lorentz factor can explain 
the formation of mini-jets with very high effec- 
tive bulk Lorentz factors. This scenario would 
avoid problems associated with very high bulk 
Lorentz factors of the jet itself, namely problems 
concerning the stati stics of detected objects 
( Henri fc Saugel [2006) and concerning tension 
between high inferred Lorentz factors on the 
order of 50 and rather slow observed pattern 
speeds in radio in terferometric observations (e.g. 
Piner et al. ( 2010| ) ) . Future work on reconnection 
might corroborate this possible link between jet 
formation the non-thermal emission. 
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